Total Cross Section, Inelasticity and Multiplicity Distributions in 

Proton — Proton Collisions. 



G. Musulmanbekov 

Multiparticle production in high energy proton - proton collisions has been anal- 
ysed in the frame of Strongly Correlated Quark Model (SCQM) of the hadron struc- 
ture elaborated by the author. It is shown that inelasticity decreases at high energies 
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and this effect together with the total cross section growth and the increasing with 
collision energy the masses of intermediate clusters result in the violation of KNO - 



. scaling. 
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(r*) ■ In inelastic hadronic interactions with multiparticle production only the fraction of colli- 

sion energy is converted into the production of secondaries. For the quantitative estimation 
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of this fraction one can use for a given interaction a characteristic, originating from cosmic 



ray physics, inelasticity, which can be defined as 
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where s - is square of cm. energy and M is the mass of intermediate system which decays 
into final produced particles. The remaining part of incident energy is carried away by 
participant's remnants - so-called leading particles. From experimental point of view more 
suitable is another definition of inelasticity 

h(s) = ^fJ^/ dy/M-^-cosliy, ( 2 ) 



where fii = w p\._ + mf is transverse mass of a produced particle of type i and drii/dy is its 
measured rapidity distribution. Fluctuation of inelasticity from event to event leads to the 
distribution P(k) with mean inelasticity < k(s) >. The energy dependence of inelasticity is a 
problem of great interest both from theoretical and experimental points of view. There is no 
consensus in the physical community on the energy dependence of < k(s) >. The decrease 
of inelasticity with energy is advocated by some authors 
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while others believe that 



inelasticity is an increasing function of energy 0, 0, Q, 0| . The question can not be answered 
by collider experiments. At ISR energies (23 -j- 60 GeV) where leading particle spectrum 
could be measured the inelasticity defined to be about 0.5. In the collider experiments at 
higher energies (SPS and Tevatron) leading particles are emitted in an extremely forward 
cone and could not be measured due to the presence of the beam pipe. Obviously, multiplicity 
distributions are connected with inelasticity distributions, and so one can study features of 
multiplicity distributions deriving the information on inelasticity or fraction of the initial 
energy converted into the particle production. As we know scaled multiplicity distributions 
exhibit KNO p scaling up to ISR energies which is violated for higher energy data (SPS) 
where they can be described approximately by the negative binomial distribution (NBD). 
And again at the most high SPS energy, 900 GeV, there is an evident deviation of multiplicity 
distributions from NBD. In this paper we demonstrate that there is a connection between 
the energetic behaviour of the shapes of scaled multiplicity distributions and inelasticity 
distributions which, in turn, relates to the effect of the total cross section growth. For this 
purpose we use geometrical considerations which are justified by the following arguments. 
First, hadrons are extended objects with the size about 1 fm and , second, at high energies 
de-Bloglie wave-length becomes small. 

Our analysis is based on the model of hadron structure so-called Strongly Correlated 
Quark Model (SCQM) which is described in Section 2. In Section 3 the model is 
applied for the calculation of proton - proton total cross sections and simulation of inelastic 
events. 



II. STRONGLY CORRELATED QUARK MODEL 

The ingredients of the model are the following. A single quark a of definite colour 
embedded in vacuum starts to polarize its surrounding that results in the formation of quark 
and gluon condensate. At the same time it experiences the pressure of vacuum because of 
zero point radiation field or vacuum fluctuations which act the quark tending to destroy 
the ordering of the condensate. Suppose that we place the corresponding antiquark in the 
vicinity of the first one. Owing to their opposite signs colour polarization fields of quark 
and antiquark interfere destructively in the overlapped space regions eliminating each other 
at most in the middle-point between the quarks. This effect leads to the decreasing of 
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condensates density in the same space region and overbalancing of the vacuum pressure 
acting on quark and antiquark from outer space regions. As a result the attractive force 
between quark and antiquark emerges and quark and antiquark start to move towards each 
other. The density of the remaining condensate around quark (antiquark) is identified with 
the hadronic matter distribution. At maximum displacement in qq— system that corresponds 
to small overlapping of polarization fields the hadronic matter distributions have maximum 
extent and magnitude. The closer they to each other, the larger is the effect of mutual 
destruction and the smaller hadronic matter distributions are around quarks and the larger 
their kinetic energies. In that way quark and antiquark start to oscillate around their 
middle-point. For such interacting qq— pair located on X axis at the distance of 2x from 
each other the total Hamiltonian is 

H = ^ _ £2)1/2 + (1 - /?2)l/2 + ^?<?( 2X )' ( 3 ) 

were m-g, m q - current masses of valence antiquark and quark, (5 — j3{x) - their velocity 
depending on displacement x and Vg q - quark-ant iquark potential energy at separation 2x. 
It can be rewritten as 



H = 



" >7J + U(x) 



(1-/5 2 )V2 



+ 



1 U(x) 



H^ + H,, (4) 



^l- 02)1/2 

were U(x) = \Vq q (2x) is potential energy of quark or antiquark. Quark (antiquark) with the 
surrounding cloud (condensate) of quark - antiquark pairs and gluons, or hadronic matter 
distribution, forms the constituent quark. It is natural to assume that the potential energy 
of quark (antiquark), U(x), corresponds to the mass Mq of the constituent quark: 

/oo /»oo /»oo 

dz' dy' dx'p(x,r') « 2Mq(x) (5) 
-oo J —oo J —oo 

where C\ is a dimensional constant and hadronic matter density distribution, p(x, r'), is 
defined as 

p(x,r') = C 2 \<f(x,r')\ = C 2 \ <p Q (x' + x, y', z') - ^q{x' - x, y', z') \ . (6) 

Here C 2 is a constant, ^>q and <p-Q are density profiles of the condensates around quark and 
antiquark located at distance 2x from each other. Here, we consider that the condensates 
around quark and antiquark have opposite colour charges. They look like compressive stress 
and tensile stress (around defects) in solids. The generalization to a three-quark system 
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in baryons is performed according to S77(3) co ior symmetry: in general, pair of quarks have 
coupled representations 

3®3 = 6©3 (7) 

in SU (3) co ior and for quarks within the same baryon only the 3 (antisymmetric) representa- 
tion occurs. Hence, the antiquark can be replaced by two correspondingly coloured quarks 
to get a colour singlet baryon and the destructive interference takes place between colour 
fields of three valence quarks (VQs). Putting aside the mass and charge differences and 
spins of valence quarks we may say that inside the baryon three quarks oscillate along the 
bisectors of the equilateral triangle. Therefore, keeping in mind that quark and antiquark 
in mesons and three quarks in baryons are strongly correlated, we can consider each of 
them separately as undergoing oscillatory motion under the potential (JHJ) in 1+1 dimension. 
Hereinafter, we consider VQ oscillating along X— axis, and Z— axis is a perpendicular to 
the plane of oscillation XY . Density profiles of condensates around VQs are taken in the 
gaussian form. It has been shown in papers Q| that the wave packet solutions of time 
dependent Schrodinger equation for a harmonic oscillator move in exactly the same way 
as corresponding classical oscillators. These solutions are called "coherent states". This 
relationship justifies our semiclassical treatment of quark dynamics. 

We specify the mass of constituent quark at maximum displacement as 

M g(g) 0w) = \ (^±^)« 360 MeV, 

where m& and are masses delta-isobar and nucleon respectively. The parameters of 
the model, namely, maximum displacement, x max , and parameters of the gaussian function, 
a x,y,z, for hadronic matter distribution around VQ are chosen inside the following corridors: 

x m ax = 0.64 -=- 0.66 fm, a XjV = 0.24 -=- 0.28 fm, a z = 0.12 + 0.20 fm. (8) 

They are estimated by comparison of calculated and experimental values of inelastic cross 
sections, 0i n (s), and the inelastic overlap function G m (.s, b) iorpp— and pp— collisions (see the 
next section). The current mass of the valence quark is taken to be 5 MeV. The behaviour 
of potential © evidently demonstrates the relationship between constituent and current 
quark states inside a hadron (Fig. QJ. At maximum displacement quark is nonrelativistic, 
constituent one (VQ surrounded by the condensate), since the influence of polarization 
fields of other quarks becomes minimal and the VQ possesses the maximal potential energy 
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corresponding to the mass of the constituent quark. At the origin of oscillation, x = 0, 
antiquark and quark in mesons and three quarks in baryons, being close to each other, 
have maximum kinetic energy and correspondingly minimum potential energy and mass: 
they are relativistic, current quarks (bare VQs). This configuration corresponds to so-called 
"asymptotic freedom" . In the intermediate region there is increasing (decreasing) of the 
constituent quark mass by dressing (undressing) of VQs due to decreasing (increasing) of 
the destructive interference effect. The evolution of colour charge density profiles of quark 
- antiquark pair during the half-period of oscillation is shown in Fig. |21 Here, we suppose 
that the quark colour charge is positive and the antiquark color charge is negative. 

The proposed dynamical picture meets local gauge invariance principle. Indeed, destruc- 
tive interference of color fields of quark and antiquark in mesons and three quarks in baryons 
depending on their displacements can be treated as a phase rotation of wave function of a 
single VQ in colour space ip c on angle 9 depending on the displacement x of the VQ in the 
coordinate space 

- e^V c (x). (9) 

The color phase rotation, in turn, leads to VQ dressing (undressing) by quark and gluon 
condensate that corresponds to the transformation of a gauge field 

A^x) ^ A^x) + d^9(x). (10) 

Here, we drop colour indices of A^(x) and consider each quark of specific colour separately 
as changing its effective colour charge, g8(x), in color fields of other quarks (antiquark) 
due to the destructive interference. Thus, gauge transformations Q, (fTUj) map internal 
(isotopic) space of a colored quark onto the coordinate space. On the other hand, this 
dynamical picture of VQ dressing (undressing) corresponds to the chiral symmetry breaking 
(restoration). Due to this mechanism of VQs oscillations the nucleon runs over the states 
corresponding to the certain terms of the infinite series of Fock space 

I B) = c x | qiq 2 q-i) + c 2 | q^qiqq) + c 3 | q^q^g)-- (H) 

The proposed model has some important consequences. Inside hadrons the valence quarks 
and accompanying them gluons and quark-antiquark pairs, as well, are strongly correlated. 
Nucleons are nonspherical object: they are flattened along the axis perpendicular to the 
plane of quarks oscillations. 
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From the form of the quark potential (Fig. one can conclude that dynamics of VQ 
corresponds to nonlinear oscillator and VQ with its surrounding can be treated as a nonlinear 
wave packet. Moreover, our quark - antiquark system turned out to be identical to so-called 
"breather" solution of (nonlinear) sine-Gordon (SG) equation^]. SG equation in (1+1) 
dimension in the reduced form for scalar function <p(x, t) is given by 



□0(x, t) + sin 0(x, t) = 0, 



(12) 



where x, t-dimensionless. Breather is a periodic solution representing a bound state of a 
soliton-antisoliton pair which oscillates around its center of mass: 



4>br{%, t) = 4 tan 
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sinh 


Ut/y/il ~ U 2 






u cosh 


x/y/(l — U 2 





(13) 



where u is 4-velocity. During the oscillations of the soliton-antisoliton pair their density 
profile 



lfi br (x,t) 



dcpbrix, t) 

dx 



(14) 



evolves like our quark-ant iquark system, i.e. at the maximal displacement the soliton and 
antisoliton are emphasized maximal and at the minimum displacement they "annihilate" 
(Fig. El). This similarity is not surprising because our quark-antiquark system was formu- 
lated in close analogy with the model of dislocation-antidislocation which in the continuous 
limit is described by the breather solution of SG equation Q]. It can be shown that soli- 
ton, antisoliton and breather obey relativistic kinematics, i.e. their energies, momenta and 
shapes are transformed according to Lorentz transformations. Since the above considera- 
tion of quarks as solitons is purely classical the important problem is to construct quantum 
states around them. Although the soliton solution of SG equation looks like an extended 
(quantum) particle the relation of classical solitons to quantum particles is not so trivial. 
The technique of quantization of the classical solitons with the usage of various methods 
has been developed by many authors. The most known of them is semiclassical method of 
quantization (WKB) which allows one to relate classical periodic orbits (breather solution 
of SG) with the quantum energy levels 14^ . 

Hereinafter, we adhere to our semiclassical model (SCQM) applying it for analysis of 
cross sections and the multiparticle production in hadron-hadron collisions. 
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III. HADRON - HADRON COLLISIONS 



Different configurations of the quark content inside a hadron realized at the instant of 
the collision result in different types of reactions. The probability of finding any quark 
configuration inside the hadron is defined by the probability of VQ's displacement in a 
proper frame of the hadron: 



where E q is a total energy of the valence quark (antiquark) and a constant A can be derived 
from the normalization condition 



Configurations with nonrelativistic constituent quarks ( x ~ a; max ) in both colliding hadrons 
lead to soft interactions with the nondiffractive multiparticle production in central and 
fragmentation regions (Fig. 0^). The hard scattering with the jet production and the 
large angle elastic scattering take place when configurations with current VQs (x ~ 0) in 
both colliding hadrons are realized (Fig. 0Jo). The near current quark configuration inside 
one of the hadrons and constituent quark configuration inside the second one result in 
single diffraction scattering (Fig. 0J:). And at last intermediate configurations inside one or 
both hadrons are responsible for semihard and double diffractive scattering (Fig. The 
same geometrical consideration can be applied to deep inelastic scattering processes if one 
assumes that a real or virtual photon converts into the vector meson according to the vector 
dominance model. 

First, we apply our model for the calculation of proton-proton and antiproton-proton 
cross sections at high energies and demonstrate that the growth of the cross section with 
energy is caused by predominantly increasing contribution of peripheral interactions that, 
in turn, leads to the decreasing of inelasticity of collisions. Then we will show that the 
energetic behaviour of inelasticity distributions governs the energetic behaviour of scaled 
multiplicity distributions. 




(15) 




(16) 
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A. Cross Sections 

To calculate cross sections we used an impact parameter representation, namely Inelastic 
Overlap Function (IOF), which can be specified via the unitarity equation 

2Imf(s,b) = \f(s,b)\ 2 + G m (s,b), (17) 

where /(s, 6)— elastic scattering amplitude and G in (s,b) is IOF. IOF is connected with in- 
elastic differential cross sections in impact parameter space: 

\da-Jdb 2 ) =G in (s,6). (18) 

7T 

Then the inelastic, elastic and total cross sections can be expressed via IOF as 

(T ta {s) = [ G in (s,b)d 2 b, (19) 



OtotlS) 



2 



1 - y/1 - G in (s, b) d 2 b, (20) 



l-y/l-G- m (s,b)\d 2 b. (21) 

Since IOF relates to the probability of inelastic interaction at a given impact parameter, we 
carried out Monte Carlo simulation of inelastic nucleon-nucleon interactions. The Inelastic 
interaction takes place at the definite impact parameter b if the produced mass meets the 
following requirement: 

M 2 F = AM PlP M TlT J P p(t) Pt (v - b) d 3 v > (M 2 F ) min , (22) 

where Mqf is the mass of the central "fireball" produced at the ovelapped region, pp and 
Pt are hadronic matter density distributions in projectile and target hadrons, Mp, Mp - 
masses, jp, 7t - gamma-factors of the colliding hadrons and (Mp F ) m j n - the minimal mass of 
a fireball that results in an inelastic event. This expression is the modification of assumption 
of Heisenberg on interaction of extended particles: we replaced in his original formula 
the pion mass squared (on the right hand side) by (M^, F ) min . In our previous papers this 
quantity corresponded to the transverse mass of the pion: m 2 ± = p\ + m 2 . Taking into 
account the energy dependence of the average momentum of produced particles and the 
increasing yield of minijets (as treated in what follows) we parametrize the minimal fireball 
mass as 

(M CF ) min = 0.3 + 0.03s 1/4 . (23) 
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Specifying the quark configurations in each colliding hadrons according to (|15|) we calcu- 
lated G in (s, b) for particular values of the impact parameter b and then according to (J2*T]) - 
total cross sections, cr to t- Fig. Elshows the results of the calculation for total cross sections for 
proton - proton and antiproton - proton collisions in a wide range of collision energies. One 
can see that the model with fixed parameters characterizing the geometrical size of hadrons 
describes the energetic behaviour of a to t rather well. The growth of the total cross section 
with energy coming from the growth of the inelastic cross section is due to the continuous 
tails of condensates (hadronic matter distributions) around VQs not compensated by the 
destructive interference effect inside each interacting particle. With rising collision energy 
the overlap of more peripheral parts of these tails make it possible to meet requirement (f2*2*j) 
and consequently results in the increasing effective size of the hadronic matter distribution 
inside nucleons and correspondingly the increasing radius of interactions. It can be seen 
from the comparison of IOFs for ISR and SPS energies which is given in Fig. |3 According 
to Eq. (jl8)l the difference Gf^ s — G[^ R exhibits the predominantly peripheral increase of the 
inelastic cross section (and thus of the total cross section, since <j e ij a to t is only about 20%) 
which is centered around 1 fm (Fig. IHb)- As noted by the authors of the paper 17] at high 
energies colliding nucleons become blacker, edger and larger ("BEL-effect"). The model 
gives the linear logarithmic energy dependence for total cross sections. At energies \fs < 30 
GeV calculated cross sections were corrected on contributions of Regge poles exchange by 
using Donnachie and Landshoff parametrization^^. An oscillatory motion of VQs appear- 
ing as an interplay between constituent and bare (current) quark configurations results in 
fluctuations of the hadronic matter distribution inside colliding nucleons. The manifesta- 
tion of these fluctuations is a variety of scattering processes, hard and soft, in particular, the 
process of single diffraction (SD). SD-events correspond to the constituent quark configu- 
ration inside one colliding hadron and (semi)bare quark configuration inside the other one. 
Our unified geometrical explanation of diffractive, nondiffr active and DIS processes could 
give an answer to the long standing question: what is pomeron? Historically the concept 
of "Pomeron" originating from simple Regge pole with the intercept ao = 1 transformed 
to a rather complicated object with relatively arbitrary features and smooth meaning. To 
produce the rising cross sections it must have the intercept such that a = 1 + e. The fact 
that the parameter e is universal, independent of particles being scattered in hadronic and 
DIS interactions, could say us that the nature of the cross section growth is the same for 
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all processes. Our interpretation of pomeron is a geometrical one. Both diffractive and 
nondiffractive particle productions emerge from the disturbance (excitation) of overlapped 
continuous vacuum polarization fields (gluon and qq condensate) around valence quarks of 
colliding hadrons followed by fragmentation process. The type of the interaction depends 
on quark configurations inside a colliding hadron occurring at the instant of the interaction 
and the value of the impact parameter. So, what we used to call "Pomeron" in t— channel 
is solely continuum states in s— channel and we claim that Pomeron is unique in elastic, 
inelastic (diffractive and nondiffractive) and DIS. 



According to our model the configurations with nonrelativistic constituent quarks ( x ~ 
Xmax) inside the both colliding hadrons lead to soft interactions with multiparticle production 
in central and fragmentation regions. The additional restriction by small impact parameters 
selects central collisions when hadronic matter distributions of colliding hadrons (quarks) 
overlap totally. In this case kinetic energies of colliding hadrons dissipate totally converting 
into the production of secondary particles that corresponds to collision inelasticity close to 
1 and very high multiplicity in comparison with the mean one. We will consider the soft 
interactions and nondiffractive multiparticle production, in particular. According to the 
KNO hypothesis the scaled multiplicity distributions, (n) P n (s), depend on the ratio of the 
number of particles to the average multiplicity z — nj (n) and they are energy independent. 
From our geometrical point of view such behaviour could be explained as a superposition 
of relatively narrow distributions corresponding to the particular impact parameters of the 
collisions. Indeed, the multiplicity distribution can be defined as 



where P[k(s)] is inelasticity distribution and P{n \ k) is the probability of the production 
of n particles at the given inelasticity, k. So, if the conditional probability P(n \ k) is 
sufficiently narrow then the shape of the distribution P n is defined by the shape of the 
inelasticity distribution P[k(s)}. The inelasticity distributions are strictly connected with the 
impact parameter distributions. KNO scaling holds (at least, approximately) if the impact 
parameter distributions and consequently inelasticity distributions are energy independent. 



B. Multiparticle Production in Hadronic Collisions 




(24) 
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As shown in the previous subsection the growth of inelastic and total cross sections with 
energy in hadronic collisions is due to the increasing of effective sizes of interacting hadrons. 
To make a quantitative analysis of energetic dependence of multiplicity distributions we 
performed, in the frame of SCQM, Monte Carlo simulation of inelastic proton - proton 
interactions selecting nondiffr active events. The process of simulation includes the following 
steps. 

1. Applying Heisenberg prescription (}2~2~)) we define the mass of the central "fireball" (CF) 
(Fig. C|) produced in the proton-proton collision at a particular impact parameter. 
Quark configurations inside each proton at the instant of the collision is specified 
randomly according to the probability (|15jl that allows one to fix energies and momenta 
of quarks inside both protons. Since the mass of CF is formed by the overlap of 
hadronic densities of individual constituent quarks (CQ) of colliding protons we know 
the energies and momenta of quarks in both remnants which we call, by convention, 
forward and backward "fireballs" (FF and BF). The notion "fireball" is applied by 
convention only because all fireballs, CF, FF and BF, can decay string - like manner 
and there is no sharp boundary between secondaries emitted from fireballs in the 
rapidity space for nondiffractive events. Then we calculate the effective masses of FF 
and BF (Fig. EJ): 



M; 



FF 



M, 



BF 




(25) 



(26) 



where E i and kj are energies and momenta of constituent quarks after the collision. 

2. We assume that each fireball breaks up, in general, into clusters. Here, the 
bremsstrahlung analogy is used, namely, at the instant of the collision a proton (elec- 
tron) loses the energy dumping fraction of its hadronic (electromagnetic) field by means 
of the emission of clusters (photons). To simulate the masses of the clusters we apply 
the result of the paper [19] for a cluster mass spectrum 



P{m d ) = (m c i/m )exp(-m d /m ), 



(27) 
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following from the statistical nature of the cluster emission. Our next assumption is 
that masses of clusters increase with the collision energy. This is dictated by the neces- 
sity to take into account such peculiarities of multiparticle production as the growth 
of rapidity distribution plateau, the increasing of transverse momenta of secondaries 
and the increasing yield of minijets. We parametrize the energy dependence of the 
average mass of clusters as 

(m d ) = 0.3 + 0.09s 1/4 . (28) 

Notice that we have chosen the same energetic dependence for the minimal fireball 
mass in Heisenberg prescription (}2*3*j) except for the value of the slope parameter. 

3. Given the positions of the centers of masses of each fireball in rapidity space and 
kinematically allowed rapidity (sub) spaces for the breaking up of each fireball into 
clusters we simulate the momenta for each generated cluster in the proper frame of 
the corresponding fireball. Bremsstrahlung mechanism of the fireball fragmentation 
corresponds to statistically independent emission of clusters with the limited transverse 
momenta. Therefore, we apply the cylindrical phase space model according to which 
the rapidity of i-ih cluster is defined as 

y* = en (29) 

where £ is random number uniformly distributed in the interval [0,1]. The allowed 
rapidity interval, Y i: is given by 

where (/i c z)i = ( m c/)? + v\n transverse mass of the cluster %. Moreover, the rapidity 
interval for fragmentation of the central fireball, Y i , is restricted by the requirement 

yCF < yFF _ yBF^ (3Q) 

where Y FF and Y BF are rapidities of forward and backward fireball respectively. The 
transverse momenta of the clusters are generated according to the distribution 



f(p\) oc exp(-bp 2 ± ). 



(31) 
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The energy of the remnant baryon in the proper frame of FF (BF) is defined by the 
summary energy of two quarks closest to each other in rapidity space, i.e. kinematic 
characteristics of the baryon are connected to those of "diquark" . 

4. Since our clusters are identified with minijets they should decay in a jet-like manner. 
One could assume that these minijets are formed by the fragmentation of the excited 
sea quark-antiquark pair. Hence, we can approximate the spectrum of cluster decay 
using data on the electron-positron annihilation provided that the cluster mass is 
identified with the center of mass energy of an electron and positron: m c \ = i/i e+e _. 
One can apply for this purpose any of appropriate Monte - Carlo generators. The axis 
of the decaying jet is generated to be directed isotropically. And, at last, the model 
meets energy - momentum conservation requirements for all products of a reaction. 

Calculated in such a manner the multiplicity distributions in the KNO form and 
the energetic dependence of the mean multiplicity for charged particles are shown in 
Figures |H1 and El Given all characteristics of produced particles in the event we can 
calculate inelasticity ki (Eq. EJ). Its distributions for pp— interactions at different 
collision energies are shown in Fig. The inelasticity distribution evolves with 

energy in such a way that its maximum position shifts to the lower values of inelas- 
ticity at higher collision energies. It means that the higher is the collision energy the 
lower is the average inelasticity (Fig. H~T| . Analyzing the multiplicity distribution at 
different energies one can see that its maximum position shifts to the lower values 
of scaled multiplicities and the contribution of high multiplicities increases while the 
collision energy increases. According to Eq. (|24|) the multiplicity distribution can 
be expressed via the conditional multiplicity distribution at particular inelasticity and 
the inelasticity distribution. The conditional multiplicity distribution at the particular 
inelasticity, in turn, is built from multiplicities of clusters emitted from forward, back- 
ward and central fireballs and multiplicities going from the clusters' fragmentation. If 
the inelasticity distribution and (average) mass of the clusters would not depend on 
the collision energy then the scaled multiplicity distributions do not depend on energy 
either and KNO - scaling takes place. The shift of the position of the maximum of 
the inelasticity distribution with the energy growth shifts the position of the scaled 
multiplicity distribution. On the other hand, masses of clusters growing with collision 
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energy lead to the narrowing of available rapidity space and, consequently, to the vi- 
olation of Feynman scaling. This effect most obviously exhibiting at the inelasticities 
close to 1 cause the lift of a multiplicity distribution tail at high multiplicities. To sum- 
marize we claim that the mean inelasticity decreases with energy, and the violation of 
KNO - scaling is a consequence of the growth of inelastic and total cross sections and 
of masses of emitted clusters with energy. 

This research was partly supported by the Russian Foundation of Basics Research, 
grant 01-07-90144. 
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FIG. 1: a) Potential energy of valence quark and mass of constituent quark; b) "Confinement" 
force. 
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FIG. 2: Evolution of color charge density profile, ip, in quark-antiquark system during half-period 
of oscillations; d = 2x - distance in fermi between quark and antiquark depicted as dots. 
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FIG. 3: Evolution of breather, <j), and its energy density profile, ip, during a half-period of oscilla- 
tion. Scales are abitrary. 




FIG. 4: Different quark configurations realized inside colliding nucleons at the instant of collision. 
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FIG. 6: Plot of overlap functions at ISR and SPS energies, (a), and the difference of overlap 
functions — G^ R , (b), as a function of impact parameter b. 
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FIG. 7: Fireball picture of multiparticle production in proton - proton collision, pi and ki are 
momenta of constituent quarks before and after collision, respectively; x\ and X2 - fractions of 
protons momenta forming the central fireball. 
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FIG. 8: The multiplicity distributions of c harg ed particles in pp— and pp— collisions. Data are 
taken from electronic data base HEPDATA 161. 
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FIG. 10: Inelasticity distributions for pp— and pp— collisions calculated according to Eq. Pj). 
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FIG. 11: The average inelasiticity plotted as a function of the energy for pp— and pp— collisions. 
Boxes are compilation of data given in paper 20]. 



